Two-dimensional dark soliton in the nonlinear Schrodinger equation 
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Two-dimensional gray solitons to the nonlinear Schrodinger equation are numerically created by 
two processes to show its robustness. One is transverse instability of a one-dimensional gray soliton, 
and another is a pair annihilation of a vortex and an antivortex. The two-dimensional dark solitons 
are anisotropic and propagate in a certain direction. The two-dimensional dark soliton is stable 
against the head-on collision. The effective mass of the two-dimensional dark soliton is evaluated 
from the motion in a harmonic potential. 
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I. INTRODUCTION 

The nonlinear Schrodinger equation is a typical soliton equation. Solitons in optical fibers are described by the 
nonlinear Schrodinger equation Recently, solitons, dark solitons and vortices were found in the Bose-Einstein 
condensates 0, 0, 0, IE • The matter- wave solitons are also described with the nonlinear Schrodinger equation 
or the Gross-Pitaevskii equation. The one-dimensional nonlinear Schrodinger equation is an integrable system, and 
the time evolution can be solved exactly. The two-dimensional nonlinear Schrodinger equation has been intensively 
studied but the behaviors are not understood completely. We study mainly a two-dimensional dark soliton to the 
two-dimensional nonlinear Schrodinger equation. The model equation is written as 
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where the Planck constant % is set to be 1, the coefficient g of the nonlinear term can be set to +1 or -1, and the 
mass parameter m is assumed to be m = 1 for the sake of simplicity. There is a two-dimensional soliton called the 
Townes soliton in the two dimensional nonlinear Schrodinger equation in the focusing case of g = —1, but the Townes 
soliton is unstable. The peak amplitude of the localized solution increases and it leads to divergence, if the norm 
of the localized solution is larger than the norm of the Townes soliton. Hereafter, we consider the defocusing case 
\ of g — 1. There is a family of one-dimensional dark soliton to the two-dimensional nonlinear Schrodinger equation, 
\Q ■ which has a form 



</>(x,y) = {iB + At&nhA{y - vt)} exp(-i/it), (2) 



where A 2 + B 2 — 1 and v = B. In the case of B = 0, the minimum value of \(j)\ becomes zero and it is called a black 
soliton. For 7^ 0, the minimum value of is \B\ and then the dark soliton is called a gray soliton. The gray 
soliton propagates with velocity B. Because the nonlinear Schrodinger equation is invariant with respect to the scale 
transformation: <\> — > a<p : x — > x/a, y — > y/a, t — > t/a 2 , qb = a{iB + Atanh Aa(y — vat)} ex.p(—ifia 2 t) is also a dark 
soliton. The velocity of the dark soliton is scaled as v — > av. There is another family of two dimensional solutions 
called vortex, which has a form qb = U(r) exp(zm#), where r is the radius from the vortex center, is the phase angle 
around the vortex center, and m is a nonzero integer. The vorticity of this vortex solution is 2i\m. The quantized 
vorticity is characteristic of the quantum fluid. The amplitude U (r) satisfies 



2 \dr 2 r dr T 2 

The amplitude U(r) at the vortex center is zero, i.e., U(0) = 0, and U(r) = 1 for r — > oo. The vortex with 
m = ±1 is called a fundamental vortex. The vortex with amplitude a is obtained by the scale transformation as 
4> = aU(ar) exp(zra0). A two-dimensional dark soliton, which has a two-dimensional hole around the center, might 
be interpreted as a kind of vortex solution corresponding to m = 0. However, it is known that there is no two- 
dimensional isotropic black soliton of the form (j)(x,y) = J7(r), which has a property of U(0) =0 and v = 0. Jones 
and Roberts found numerically two-dimensional special solutions called the rarefaction pulses, which correspond to 
two-dimensional gray solitons @, Ej- In this paper, we study a few natural creation processes of the two dimensional 
gray solitons. Then, we show the stability of the two-dimensional dark soliton against the head-on collision. Finally 
we evaluate the effective mass of the two-dimensional dark soliton from the motion in a harmonic potential. 
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FIG. 1: (a) Time evolution of \</>(x,ymin) \ for Oo = 2tv ■ 45/360. (b) Time evolution of vortex centers of the vortex and antivortex 
pair for O = 2tt -45/360. 




FIG. 2: Two dimensional dark soliton for Oo = 2tt/9. (a) Time evolution of \(j>(x,ymin)\> (b) Three-dimensional plot of the 
two-dimensional dark soliton. (c) Contour plot of the two-dimensional dark soliton at t = 56 starting from the initial condition 
with Oo = 2tt/9. 



II. INSTABILITY OF A ONE-DIMENSIONAL DARK SOLITON AND CREATION OF A 

TWO-DIMENSIONAL DARK SOLITON 

The one-dimensional dark soliton is linearly unstable for the modulation in the ^-direction. It is observed in 
numerical simulations and experiments that the one-dimensional dark soliton breaks up into vortex- antivortex pairs 
due to the development of the transverse instability 0, 0, . We have performed numerical simulations again by 
changing the B value of the initial dark soliton. The system size is set to be L x xL y = 25 x 200. We performed numerical 
simulations with the spit-step Fourier method of 256 x 2048 modes. The periodic boundary conditions are imposed. 
The initial conditions are assumed to <fi = iB — 6cos{27r(x — L x /2)/L x }sech{A(y — L y /8)} + AtdLiih{A(y — L y /8)} 
for y < Ly/2 and = iB - 6cos{27r(x - L x /2) / L x }sech{A(7L y /8 - y)} + A ta,nh{A(7L y /S - y)} for y > L y /2, 
where A = sin#o, B = cos#o, and b = 0.15 is used. That is, two one-dimensional dark solitons are initially set at 
y = L y /8 and 7/8L y , and transverse modulations are overlapped on the one-dimensional dark soliton solution. The 
black soliton corresponds to the case of #o = tt/2. As is decreased, the minimum value of the dark soliton becomes 
larger and the propagation velocity increases. The initial condition is mirror-symmetric with respect to y = L y /2 and 
x = L x /2. For Oo > 2tt ■ 42/360, the transverse modulation increases and a pair of vortex and antivortex are created. 
Figure 1(a) displays a time evolution of \</>(x , y m i n )\ for 0o = 2tt ■ 45/360. Here, ymin is the y coordinate where \<p\ 
takes the minimum value. The transverse modulation grows, and two minimum points appear, which corresponds to 
the appearance of a vortex pair. Figure 1(b) displays a time evolution of the x-coordinate of the vortex centers for 
#o = 27r • 45/360. The positions of the vortex centers were determined by the phase singularity points of (f>(x, y) in our 
numerical simulations. The vortex pair is created at x = L x /2 = 12.5, and the vortex pair is rapidly separated away. 

For 0o < 2tt - 41/360, the growth of the transverse modulation saturates in time and a pair of vortex is not created 
as shown in Fig. 2(a). The transverse modulation leads to the formation of a two-dimensional dark soliton. Figure 
2(b) is a three dimensional plot of |</>| at t = 56 for #o = 27r- 40/360. The minimum value of in the two-dimensional 
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FIG. 3: Minimum value B of \<j>\ as a function of 0o(deg.). (b) Velocity in the y direction as a function of 0o(deg/ 



dark soliton is about 0.08. We defines Bq as the minimum value of in the two-dimensionl dark soliton, which 
is different from the minimum value B of \(j)\ in the initial one-dimensional dark soliton. The nozero value of Bq 
implies a gray soliton. It is propagating in the ^/-direction. Figure 2(c) displays a contour map of \(f)\ at t = 56. The 
two-dimensional dark soliton is anisotropic. The width in the y direction is shorter than that in the x direction. 

We have calculated the minimum value Bo and the velocity v y in the ^/-direction as a function of the parameter #o 
in the initial condition. The results are shown in Fig. 3. The minimum value Bq increases almost linearly from 0, 
as #o is decreased from the critical value of the formation of the two-dimensional dark soliton. The velocity v y also 
increases almost linearly as Oo is decreased. The velocity v y is approximately expressed as v y ~ 0.38 • Bo + 0.59. This 
is a relation comparable to the relation v y = B for the one-dimensional dark soliton. There is no two-dimensional 
dark soliton with velocity v y < 0.59. 



III. MERGING OF A PAIR OF VORTICES AND CREATION OF A TWO-DIMENSIONAL DARK 

SOLITON 

The two-dimensional dark soliton might be also interpreted as a merged state of the vortex and antivortex. We 
have investigated time evolutions of vortex and antivortex pairs. The vortex pair interacts with each other. If the 
distance of the two vortices is sufficiently large, the vortices are expected to move as point vortices. The equation of 
motion of the ith point vortex among N vortices is generally described as [l3| 

dt ~ dY % " dt ~ dX % " l ~ 2tt ' W 

where Xi and Y{ are the x and y coordinates of the vortex center of the ith vortex, Yj is the vorticity of the jth 
vortex and r^- = yj {Xi — Xj) 2 + (Yl — Yj) 2 . In the simplest case of a single vortex- antivortex pair with quantized 
vorticity Yj = ±2tt, the equation of motion of the vortex centers obeys dYo/dt = l/(2xo) and dXo/dt = 0, if a 
vortex and an antivortex are initially set at (L x /2 — xo,L y /2) and (L x /2 + xo,L y /2). The vortex pair propagates 
in the ^/-direction, but does not move in the x-direction. However, if the separation of the two vortices is sufficiently 
small, the assumption of the point vortex is not good. The typical length scale is about 1/a for the vortex solution 
<j) = aU(ar) exp(±z#). The assumption of the point vortex becomes worse, if the distance between the two vortices 
is smaller than 2/ a. We have performed numerical simulations of the motion of a vortex- antivortex pair. A vortex 
and an antivortex with a = 0.347 are used as an initial condition. We have confirmed that the vortex pair propagates 
in the y direction with nearly v y = 1/(2xq) for large xq. However, if the initial distance of the vortex pair is smaller 
than xo c = 2.3, the vortex pair approaches each other and it leads to the pair annihilation. We have also observed 
complicated behaviors such that the vortex pair reappear in a certain time after the pair annihilation. 

We investigate the interaction among four vortices more in detail. If two vortices are set at (— Xq, Yq) and (Xq, —Yq) 
and two ant ivort ices are set at (Xq,Yq) and (— Xo, — Yb), the four vortices interact with one another. If the mirror 
symmetries with respect to x — and y = are satisfied during the time evolution, the equation of motion of the 
vortex center (Xq, Yq) of the vortex in x > and y > is expressed as 



dY 1 X 

(5) 
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FIG. 4: Time evolutions of (a) x— coordinate X and (b) y— coordinate Y of vortex centers, 
(a) (b) 
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FIG. 5: (a) Time evolution of \</>(L x /2,y)\ for x = 4.2. (b) Contour plot of \</>(x,y)\ at t = 200 for x = 4.2. 



from the theory of point vortex [Tj|- The solution of the point vortices satisfies the relation 4XqF 2 = a 2 (Xg + F 2 ), 
where a is an arbitrary constant. If the initial conditions are Xq(0) = Yq(0) = xo, a is equal to -\/2xo, and then Xo(t) 
decreases to X0/V2 from Xo(0) = #o, and Yo(t) increases to infinity from Yb(0) = xo, as t — > 00. The distance between 
the vortex and the antivortex approaches a smaller value V2xo owing to the interaction among the four vortices, and 
the final velocity of Yq takes a faster value of v y = 1 / (v^^o) ? compared to a single vortex pair. We have performed 
direct numerical simulations of the two pairs of vortex and antivortex. Four vortices are set at (L x /2 — xo,xo), 
(L x /2 + xo, ^o), (L x /2 — xo, L y — xo), and (L x /2 + xo, L y — xq) as an initial condition. The two pairs of vortex and 
antivortex are located near the boundaries y = and y = L y , and they are set to be mutually mirror-symmetric with 
respect to the lines y = and y = L y under the periodic boundary conditions. The four vortices therefore interact 
with one another through the boundaries at y = and y = L y . The amplitude of the vortices is a = 0.347, and the 
system size is [L x x L y ) = (100,200). Figure 4(a) displays time evolutions of the x coordinates of the vortex and 
antivortex in the region of y < L y /2 for three initial values xq — 4.2, 4.9 and 5.95. For xo = 4.2, the vortex- antivortex 
pair is annihilated at t = 60.1. The reappearance of the vortex pair was not observed. The critical value xo c — 4.3 
of xo for the pair annihilation is larger than the critical value x$ c = 2.3 in the case of a single vortex pair. This is 
owing to the interaction among the four vortices. For xo = 4.9, the vortex pair approaches once, but they separate 
away. For xq = 5.95, the vortex pair approaches initially but the distance keeps about 7.5 for large t, which is close 
to the value \[2x§ ~ 8.4 expected from the motion of the point vortex. Figure 4(b) displays time evolutions of the 
y coordinates Yq of the vortex center for the same initial values xo — 4.2, 4.9 and 5.95. For xq = 5.95, the vortex 
pair moves with a constant velocity v y ~ 0.135, which is close to the velocity 1/7.5 = 0.133 of a point vortex pair. 
As xo is decreased, the distance between the vortex pair is decreased, and therefore, the velocity in the y direction is 
increased. For x = 4.9, the average distance between the vortex pair for 150 < t < 200 is 4.4, and the velocity v y is 
about 0.195, which is slightly smaller than 1/4.4 = 0.22. The point vortex approximation might become worse, if the 
distance between the vortex pair is small. For xo = 4.2, the vortex pair is annihilated at t ~ 60.1, and the curve of 
Yo stops at the time. 

Figure 5(a) displays a time evolution of \(/>(L x /2,y)\ for xq — 4.2. The merging of the vortex pair leads to the 
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FIG. 6: Head-on collision of couter-propagating two-dimensional dark solitons starting from one-dimensional dark solitons with 
0o = 2tv/9. Time evolution of \cj)\ at the cross section x — L x /2 is shown. 



formation of a two-dimensional dark soliton. The two-dimensional dark soliton propagates with a constant velocity 
v y ~ 0.26. Figure 5(b) displays a contour plot of \<j)(x, y)\ at t = 200 for xo = 4.2. An anisotropic structure of the dark 
soliton is seen also in this case. The minimum value of \(j)\ of the two-dimensional dark soliton is nearly 0.14. The 
minimum value corresponds to Bo = 0.14 • 1/0.347 = 0.403 for the vortex with a = 1. The velocity of the dark soliton 
with a = 1 and Bq = 0.403 is estimated as 0.743 from the relation v y ~ 0.38 • Bq + 0.59 obtained from Fig. 3. The 
corresponding velocity of the dark soliton for a = 0.347 is expected to be v y = 0.743a = 0.743 • 0.347 ~ 0.26 from the 
scaling relation, which is nearly the same value as the velocity measured in the direct numerical simulation. It implies 
that the two-dimensional dark soliton obtained by the annihilation of a vortex pair is the same as the two-dimensional 
dark soliton obtained by the transverse instability of a one-dimensional dark soliton. 



IV. COLLISION OF TWO-DIMENSIONAL DARK SOLITONS 



We study the head-on collision of two two-dimensional dark solitons, to show the stability of the two-dimensional 
dark soliton. The system size is again L x x L y = 25 x 200. As an initial condition, two one-dimensional dark solitons 
of #o = 2tt/9 are set at y = L y /S and 7/8L y , which are mirror-symmetric with respect to y = L y /2 = 100. The 
time evolution of <fi(L/2,y) of two dark solitons is shown in Fig. 6. The one-dimensional dark solitons evolve into 
two-dimensional dark solitons after a while. The two-dimensional dark soliton in y < L y /2 = 100 propagates in the 
+y direction and it collides with its mirror-symmetric two-dimensional dark soliton propagating in the —y direction 
at t ~ 120. The two dark solitons pass through each other after the collision. It implies that the two-dimensional 
dark solitons are stable against the head-on collision. That is, the two-dimensional dark soliton has a property of 
"soliton". 



V. MOTION OF A TWO-DIMENSIONAL DARK SOLITON IN A HARMONIC POTENTIAL 

We study a dynamical behavior of the two-dimensional dark soliton in a harmotic potential. It is known that the 
one-dimensional dark soliton behaves like a classical- mechanical particle in a weak harmonic potential U(y) = ky 2 /2, 
and the center y of the dark soliton obeys the equation 

d2 y t 

However, the effective mass ra e // is not m = 1, which is simply expected from Eq. (1), but nearly 2 0, 0, Hf^ . 
On the other hand, the center of the bright soliton and the vortex obey the equation of motion with m e ff = 1 in a 
weak harmonic potential The motion of the background density needs to be taken into considered to understand 
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FIG. 7: (a) Time evolution of \<f>(L x /2,y)\ of a two-dimensional dark soliton in a harmotic potential U(y) = 0.00025 • y 2 . (b) 
Time evolution of the minimum value Bo. (c) Time evolution of the x coordinates Xq of a vortex pair, (d) Velocity v y of the 
two-dimensional dark soliton in the y direction. The dashed line is v y = 0.78 cos(0.018£). 

the effective mass different from 1. The difference of ra e // and m = 1 is also related to the fact that the B- value 
and therefore the velocity v y = B can take an arbitrary value in the dark soliton solution (2). The £?o-value in 
the two-dimensional dark soliton also takes an arbitrary value, and therefore the velocity v y changes continuously 
according to the Bq- value as shown in Fig. 3. We can therefore expect that the effective mass of the two-dimensional 
dark soliton might be different from 1. To check the possibility, we have studied a model equation: 

*^ = -Q) V 2 0+|0| 2 + t/(2/)0, (6) 

where U(y) = k/2(y-L y /4) 2 for y < L y /2 and U(y) = k/2(y - 3L y /4) 2 for y > L y /2 with k = 0.0005. The potential 
also satisfies the mirror-symmetry with respect to y — L y /2 = 100. As an initial condition, we have set a two- 
dimensional dark soliton at y = L y /4 = 50, which was obtained due to the transverse instability of a one-dimensional 
dark soliton of 6o = 27r/360 • 34. Figure 7(a) displays a time evolution of \<fr\ at the cross section of x = L x /2. 
Figure 7(b) displays the time evolution of the minimum value Bq of \(j>\. The £?o-value decreases with t owing to 
the harmonic potential. When the Bq value reaches 0, the two-dimensional dark soliton cannot exist and a pair of 
vortices is created from the two-dimensional dark soliton. Figure 7(c) displays a time evolution of the x coordinate 
of the vortices. The creation of the vortex pair occurs at t ~ 53. Figure 7(d) displays the time evolution of the 
velocity v y of the two-dimensional dark soliton in the ^-direction before the creation of the vortex pair. The velocity 
v y is approximated as v y ~ 0.78 cos(0.0182£). For a classical-mechanical particle with mass m e ff in the harmonic 
potential, the velocity obeys v y — v§ cos{(/c/m e //) 1//2 £}. Figure 7(d) implies that the two-dimensional dark soliton 
behaves like a classical- mechanical particle in the harmonic potential, however, the effective mass m e ff is evaluated 
as m e ff ~ /c/0.0182 2 ~ 1.5. Thus, we have confirmed that the effective mass of the two-dimensional dark soliton is 
different from 1, although we have not yet succeeded in calculating the effective mass theoretically. 

VI. SUMMARY 

We have numerically created two-dimensional dark solitons using two kinds of processes. One is the transverse 
instability of one-dimensional dark solitons, and the other is the pair annihilation of a vortex and an antivortex. The 
two-dimensional dark solitons are propagating gray solitons. The dark soliton is anisotropic in space , i.e. the width of 
the amplitude is smaller in the direction of the propagation. We have shown that the two-dimensional dark soliton is 
stable against the collision. Finally, we have investigated the dynamical behavior of the two-dimensional dark soliton 
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in a weak harmonic potential and found that the effective mass is different from 1. 
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